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Abstract 

We have measured the branching fractions for the exclusive decay modes B 
Xci(2)-f^(-^*) using a 140 fb~^ data sample collected by the Belle detector at the 
KEKB asymmetric-energy e+e^ collider. The measured branching fractions for 
B+ XciK^, 50 ^ XciK°, XciK*° and B+ XciK*+ decay modes 

are (4.5 ± 0.2 ± 0.5) x 10"^, (3.5 ± 0.3 ± 0.5) x 10"^, (3.1 ± 0.3 ± 0.7) x 10"^ and 
(4.1 lb 0.6 lb 0.9) X 10~^, respectively, where the first error is statistical and the 
second error is systematic. We do not observe statistically significant signals for the 
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B Xc2K{K*) decay modes and set upper limits at the 90% confidence level. We 
also study the helicity distribution for B XciK* decay mode and show that the 
longitudinal polarization component is dominant. 

Keywords: S-meson, Charmonium, Branching Fractions, Polarization 

PACS: 13.25.Hw, ll.30.Er 



1 Introduction 



Decays of B mesons to two-body final states including a charmonium meson 
have played a crucial role in the observation of CP- violation in the B meson 
system [1]. These decay modes also provide a sensitive laboratory for studying 
strong interaction (QCD) effects in a heavy meson system. One widely used 
approximation to handle QCD effects in heavy meson decays is "factoriza- 
tion", where it is assumed that participating quarks form hadrons without 
subsequent transfer of quantum numbers between the hadrons. In the fac- 
torization limit, two-body B decays of the type B — > XcoX and Xc2X are 
suppressed by angular momentum and vector current conservation [2,3]. Belle 
and BaBar have reported the observation of exclusive B^ — > XcqK^ [4,5,6] 
and the measurement of inclusive B — > Xc2X decays [7]. A large B — > XcqK 
signal yield can be explained by either factorization violation or rescattering 
of D'f^D^*^ in the final state [8] ^ . This mechanism also predicts a branching 
fraction for B Xc2K comparable to that for B XcoK. The measure- 
ment of new branching fractions of both factorization allowed B — > XciK* 
and factorization inhibited B — > Xc2K* decays provides valuable information 
for understanding of Xc production in B meson decay. In this paper, we report 
a study of exclusive decays B — > Xci(2)K{K*) for both charged and neutral B 
mesons, where both Xci and Xc2 are reconstructed from the same final state. 

These measurements are based on a 140 fb~^ data sample which contains 
152 million BB pairs, collected with the Belle detector [9] at the KEKB 
asymmetric-energy e+e" (3.5 GeV on 8 GeV) collider [10] operating at the 
T(45') resonance. 



on leave from Nova Gorica Polytechnic, Nova Gorica, Slovenia 
^ Throughout this paper, when a mode is quoted the inclusion of charge conjugate 
mode is implied. 
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2 Belle detector 



The Belle detector is a large-solid-angle magnetic spectrometer. Closest to the 
interaction point is a three- layer silicon vertex detector (SVD), followed by a 
50-layer central drift chamber (CDC), an array of aerogel threshold Cerenkov 
counters (ACC), a barrel-like arrangement of time-of-fiight (TOF) scintilla- 
tion counters, and an electromagnetic calorimeter (ECL) comprised of Csl(T/) 
crystals. These subdetectors are located inside a superconducting solenoid coil 
that provides a 1.5 T magnetic field. An iron flux-return located outside of the 
coil is instrumented to detect K^^ mesons and to identify muons. The detector 
is described in detail elsewhere [9]. 



3 Event selection 

Events with 5-meson candidates are first selected by applying general hadronic 
event selection criteria. These include a requirement on charged tracks (at 
least three of them should originate from an event vertex consistent with the 
interaction region), a requirement on the reconstructed center-of-mass (CM) 
energy 

(^CAf > o.2v^, where ^/s is the total CM energy), a requirement on 
the longitudinal (;s-direction) component of the reconstructed CM momentum 
with respect to the beam direction ([p^ '^^| < CSy^/c) and a requirement on 
the total ECL energy (O.ly^ < -5^^^ < 0.8i/i) with at least two energy 
clusters. To suppress continuum background, we also require that the ratio 
of the second and zeroth Fox-Wolfram moments [11] be less than 0.5. To re- 
move charged particle tracks that are poorly measured or do not come from 
the interaction region, we require dz < 5 cm for all tracks other than those 
identified as decay daughters of Kg, where dz is the track's closest approach 
to the interaction point along the beam direction. 



4 Xci(2) meson candidates 

The Xci(2) meson states are reconstructed via the decay mode Xci(2) iJ/"^ 
with J/ip ^ where I is muon or electron. For muon tracks, identifi- 

cation is based on track penetration depth and the hit pattern in the KLM 
system [12]. Electron tracks are identified by a combination of dE/dx from 
the CDC, E/p (E is the energy deposited in the ECL and p is the momentum 
measured by the SVD and the CDC) and shower shape in the ECL [13]. In or- 
der to recover di-electron events in which one or both electrons have radiated 
a photon, the four-momenta of all photons within 0.05 radian of the e"*" or e~ 
directions are included in the invariant mass calculation. The invariant mass 
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windows arc -0.06 (-0.15) GcV/c^ < Me+e--Mj/^ < 0.036 GeV/c^ to select 
J/ip candidates in the fi^ fi^ (e^e^) channels; these intervals are asymmetric 
in order to include part of the radiative tails. Vertex- and mass- constrained 
kinematic fits are then performed for selected J/ ip candidates to improve the 
momentum resolution. The Xc states are reconstructed by combining a. J/ip 
candidate with momentum below 2 GeV/c in the CM frame with all photons 
of energy greater than 0.060 GeV in the laboratory frame; to suppress 7r° we 
veto photons that, when combined with another photon in the event, satisfy 
0.110 GeV/c^ < m^^ < 0.150 GeV/c^. The Xci{2) candidates are selected by 
requiring the mass difference (M^+^-^ — M^+i-) to he between 0.385 GeV/c^ 
and 0.431 GeV/c^ for Xci and between 0.436 GeV/c^ and 0.482 GeV/c^ for 
Xc2- The mass difference resolutions are 9.6 MeV/c^ and 12.2 MeV/c^ at Xci 
and Xc2, respectively. Mass-constrained fits are applied to all selected Xci and 
Xc2 candidates. The {Mi+£-^ — M£+£~) mass distribution is shown in Fig. 1 
together with fitted curve [7] . 



5 K or K* candidates 



The K*^ and K*~^ candidates are reconstructed from four decay modes: K*^ 
K~^7r^ or KgiT^, and K*^ K^tt'^ or Kgir^. The kaon likelihood ratio for 
a track is defined as C{K/it) = C{K)/{C{K) + C((k)). Charged kaons are 
identified by requiring C{K/'k) to be greater than 0.6. The kaon and pion 
likelihoods {C{K) and /^(vr)) are obtained by combining measurements from 
the TOF and the CDC (dE/dx) with hit information from the ACC. Tracks 
which are not identified as either a kaon or a lepton are treated as charged 
pion candidates. In the reconstruction of the decay mode Xcii2)K^, 
all well-measured tracks other than the ones from the J/ip are considered to 
be kaon candidates. In order to eliminate the systematic error from parti- 
cle identification and to retain high efficiency, no particle identification cut for 
kaons is applied. The Kg candidates are reconstructed from pairs of oppositely 
charged pions. The tt+tt^ pair is required to be displaced from the interaction 
point (IP) by a minimum transverse distance of 0.22 cm for high momentum 
(> 1.5 GeV/c) candidates and 0.08 cm for those with momentum less than 
1.5 GeV/c. The direction of the pion pair momentum must agree with the 
direction defined by the IP and the vertex displacement within 0.03 rad for 
high momentum candidates (> 1.5 GeV/c), within 0.1 rad for candidates hav- 
ing momentum between 1.5 GeV/c to 0.5 GeV/c and within 0.3 rad for the 
remaining low momentum candidates. The invariant mass of the 7r"'"7r~ pair 
candidate is required to satisfy 0.482 GeV/c^ < M^+^- < 0.514 GeV/c^. 
The vr*^ candidates are reconstructed from pairs of photons with energies 
greater than 0.060 GeV. The tt" candidates are required to have invariant 
mass 0.120 GeV/c^ < M^^ < 0.150 GeV/c^. A mass-constrained fit is per- 
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formed to obtain the momentum of the 7r°. To reduce the background from 
slow pions, we require the 7r° momentum to be greater than 0.2 GeV/c in the 
CM frame. 



6 B meson reconstruction 

We reconstruct B mesons by combining a Xci or Xc2 candidate with a charged 
or neutral K{K*) candidate. We form two independent kinematic variables: 
the beam constrained mass Mbc = ^(-Ebeam/c^)^ — ((Pxci(2) + PW(_ft:*))/c)^, and 
AE = (-E;^^i(2) + Ek{k*)) - -Ebeam, where -Bbeam is the beam energy, and Pxci(2)' 
^XcU2) ' Pk{k*) and Ek{k*) are momenta and energies in the CM frame of the 
Xc states and the K{K*) mesons, respectively. The S-meson signal window is 
defined as 5.27 GeV/c^ < Afbc ^ 5.29 GeV/c^ for all channels and the ranges 
for AE are ± 0.025 GeV for B Xci(2)^(^s), and ± 0.040 GeV for B 
Xci{2)K*{K+7r- or K'Itt+). A larger range (-0.050 GeV < < 0.030 GeV) 
is used for B — > Xci{2)K*'^{K'^7r^ or K^n^) to accommodate the wider AE 
distribution that results from energy leakage in the calorimeter for 7r° modes. 



7 Signal yield and background estimation 

To study various sources of background, we use 5.5 million B ^ ip inclusive 
Monte Carlo events generated using the event generator EvtGen [14] and a 
GEANT3-based program [15] to simulate the Belle detector response. This 
Monte Carlo sample contains decays of the type B ^ k, where /c is a three- 
body final state, e.g. XciKn or J/ipKir etc., or two-body decays, e.g. J/ipK or 
J/t/jK* or Xci{2)K{K*) etc. Based on this Monte Carlo study, the backgrounds 
for B Xci{2)K{Kg) decay modes can be categorized as J/ijj inclusive and 
feed-across backgrounds. The J/ip inclusive background occurs due to a true 
J/ip candidate accidentally combining with an unrelated 7 to fake a Xci(2)- 
Feed-across background is due to one mode being reconstructed as a different 
mode because of the similarity in event kinematics. For B Xci{2)K* decay 
modes, there is an additional background due to non-resonant decays of B 
mesons, i.e. B — > XcI(2)Ktt. Non-resonant background peaks at the same place 
as signal in the Mbc and AE distributions, but does not peak in the /^vr-mass 
distribution. We thus extract the number of signal events from a fit to the 
Mbc distribution for B — > Xci{2)K{K^) decay modes but use the i^Tr-mass 
distribution for B —>■ Xcii2)K* decay modes. Each component of background 
has been estimated separately as discussed below. 

The J/t/j inclusive background is modeled by the sidebands in the (Me+e-j — 
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Mi+i~) mass distribution of data and scaled according to the fitted area of 
the background under the Xci{2) signal region. The sidebands are defined as 
0.2 GeV/c2 to 0.3 GeV/c^ and 0.55 GeV/c^ to 0.60 GeV/c^ in the {M^+e-^ - 
Mi+i-) mass distribution. This model includes a contribution from the combi- 
natorial t^l~ background under the J/ -0 peak ( "fake J/ -0 candidates" ) , which 
is not present in the ip inclusive Monte Carlo sample; this component forms 
about 5% of the J /ilj inchisive background. As a systematic check, the ip inclu- 
sive Monte Carlo has also been used to model the J/ if] inclusive background 
and the contribution of the J/-^ inclusive background decays have been ob- 
tained by excluding the feed-across and non-resonant background decays in 
the inclusive Monte Carlo sample. A small difference in the fit between the 
two models has been included in the systematic error. 

For each decay mode i — 1, 12, the feed-across background due to other 

Xc\{2)K{K*) decay modes j — 1, 12 with j ^ i is estimated by taking the 

appropriate Monte Carlo histograms and normalizing them according to the 
signal yields observed in the modes j. The total feed-across background for 
mode i is then fixed in the fit. The fitting procedure for all twelve XciC2)K{K*) 
modes is repeated, using the signal yields from the previous iteration as input, 
until we find a stable result. 

The non-resonant background shape is obtained from the ip inclusive Monte 
Carlo sample. Due to low statistics and similar shapes of the J/ip inclusive and 
non-resonant background, it is not possible to fioat simultaneously the nor- 
mahzations of both backgrounds in the fit to data. We have little knowledge 
about the amount of non-resonant background, so when fitting the /TTT-mass 
distribution the normalization of the J / ip inclusive background is fixed, while 
it is floated for non-resonant background. The contributions from non-resonant 
background decays have been obtained from the inclusive Monte Carlo sam- 
ple by excluding the feed-across and J/ip inclusive background decays. 

The Mbc distributions for B Xci{2)K{Kg) decay modes are fitted with a 
Gaussian to describe the signal region. The width and mean value of the 
Gaussian are fixed in the fit using Monte Carlo values that are calibrated by 
the S+ — > XciK'^ decay mode in data because this mode has the highest 
statistics among the decay modes considered. The background shapes for the 
J/ip inclusive and feed- across backgrounds have been modeled as discussed 
above. The normalization parameter for the feed-across background histogram 
is fixed while it is fioated for the J/ip inclusive background histogram. Figure 
2 displays the fit to the M^c distributions for the (a) S+ — > XciK~^, (b) 
B+ Xc2K+, (c) B° XciK°s, and (d) S° Xc2K^s decay modes. The 
signal yields extracted from these distributions are listed in Table 1. For the 
B Xc2K{Kg) mode, no statistically significant signal is seen. As a check, 
the signal yields have also been obtained by fitting the AE distributions for 
the B — > Xci(2)K{Kg) decay modes and are found to be consistent with the 
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signal yields obtained from the fits to the Mbc distributions. 

A reiativistic P-wave Breit-Wigner function [16] is used to fit the K*(892) sig- 
nal region in the i^TT-mass distributions for the B — > Xci(2)K* decay modes. For 
K*(892), the width and peak position of the P-wave Brcit-Wigncr function are 
found to be consistent between data and Monte Carlo in the B J/ipK*{892) 
decay modes [17]; therefore, those parameters are fixed to the PDG values in 
the fit [18]. The shapes of the J/ip inclusive, feed-across and non-resonant back- 
grounds are modeled as explained previously. There may be some structure 
in the Xgll^SO) region in the Kn-mass distributions for the B Xci(2)K* 
decay modes. To accommodate this possibility, we have included a reiativistic 
D-wave Breit-Wigner function in the fits to the fTTT-mass distributions. The 
mean and width of the D-wave Breit-Wigner function are fixed to the nominal 
PDG values. Figure 3 shows the fits to the i^7r-mass distribution for the (a) 
50 _ .y.^^K*\K^7,-), (b) B^ ^ Xcii^*°(i^^vrO), 5+ ^ XciK*+{K+7:^), 
and (d) B^ XciK*^{K^t^^) decay modes. The fits to the Xyr-mass distri- 
butions for (a) 5° ^ Xc2K*\K+7i-), (b) 5° ^ Xc2i^*°(^svr°), (c) B+ 
Xc2K*^{K^Ti'^), and (d) B^ — > Xc2K*^ {Kg7i~^) decay modes are presented in 
Fig. 4. The signal yields extracted from these distributions are hsted in Table 
1. For the S° XciK*'^iKs'^'^) and B — »■ Xc2K* decay modes, no statistically 
significant signals are seen. The signal yields have also been checked by fitting 
the Mbc and AE distributions for the B Xci{2)K* decay modes. For the 
Mbc and AE distributions a K* is identified if the absolute difference between 
the invariant mass of an identified Ktt pair and nominal K* mass is less than 
75 MeV/c^. The signal yields obtained from the fits to the Mbc and AE dis- 
tributions arc found to be consistent with the signal yields obtained from the 
fits to the Kir-msiss distributions. 



8 Helicity study for B — > Xci^* decays 

A P VV decay may include three possible polarization amplitudes: two 
transverse and one longitudinal. The contributions of these amplitudes can 
be obtained from the helicity distribution of the K* (cos 9k*) where 9k* 
is defined as the angle between the K* momentum and the direction op- 
posite to the B momentum in the K* rest frame. The transverse polariza- 
tion amplitudes have a sin^ 9k* dependence while the longitudinal amplitude 
has a cos^ 9k* dependence. We have analyzed the helicity distributions for 
different i^vr-mass regions. The efficiency-corrected cos 9k* distribution for 
the signal region (0.821 GeV/c^ < Mk^ < 0.971 GeV/c^) is fitted with the 
function f {cos 9k*) = a, + 3b cos'^ 9k* + 2cV5ba cos 9k*- The fitted value of c 
(0.01 ± 0.10) is consistent with zero, which shows that the interference with 
the scalar amplitude is negligible. The helicity distribution for the background 
obtained from sideband region (Mktt < 0.8 GeV/c^) is flat except for a peak 
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in the last bin due to the J/%1) inclusive background (Fig. 5(a)); this distri- 
bution is significantly difi^erent from the hehcity distribution for the signal 
region as shown in Fig. 5(b). We perform a fit to the efficiency-corrected 
cosOk* distributions for the signal and sideband regions with the function 
F{cos9k*) ^ A-\- SB cos^ 9k* and obtain Asig{Asb) and Bsig{Bsh) for the K* 
signal (sideband) region. With these parameters, the background subtracted 
net K* helicity angle distribution is extracted with corresponding parame- 
ters A* and B* by estimating A* = A^ig - SAg^ and B* = Bgig - SBgy^, 
where S is the scale factor to subtract non-resonant Kn contributions and 
other backgrounds. Thus, based on A* and B*, we obtain the longitudinal 
polarization of K* . Using the longitudinal polarization amplitude Ah and the 
total amplitude Atot, the longitudinal polarization parameter is expressed as 

\A p B* I ^* 

7] = jj^^ which is also expressed in terms of A* and B*; rj = . It is 

also observed that the difference in fit with or without inclusion of the last 
peaking bin of the cos 6k* distribution for the background region is small and 
the error assigned to rj takes into account this effect. To estimate the sys- 
tematics due to the background shape uncertainty, we obtain the value of rj 
by using the upper sideband region (1.1 GeV/c^ < Mktt < 1-3 GeV/c^) for 
a background fit. The difference in values obtained using the two different 
sideband regions is found to be 4.6%, which is treated as a systematic error. 
To estimate the influence of a possible interference term (cos 9k* ) on the po- 
larization parameter value, we perform fits taking the signal events from left 
and right sides of the K* peak i.e. (0.821 GeV/c^ < Mk^ < 0.896 GeV/c^) 
and (0.896 GcV/c^ < Mktt < 0.971 GeV/c^). In spite of the observed nonzero 
values of the scalar interference term in these distributions, the change in 
longitudinal polarization is found to be within 6.8%, which is taken as a sys- 
tematic error. These two errors have been combined in quadrature to give a 
total systematic error of 8.2% on polarization parameter (77). Hence the final 
value of longitudinal polarization parameter rj from this analysis is obtained as 
Tj = 0.87 ± 0.09 ± 0.07. From this value of rj we conclude that the B — > Xci^* 
decay is dominated by longitudinal polarization with a lower limit for the lon- 
gitudinal polarization of 0.69 at the 90% confidence level. This mode has a 
larger longitudinal polarization than B — > J/t/jK* [17,19]. 



9 Measurement of branching fractions 

The branching fraction for each B — > Xci{2)K{K*) decay mode is estimated by 
dividing the observed signal yield by the reconstruction efficiency, the number 
of BB events in the data sample and the daughter branching fractions. We 
have used the daughter branching fractions published in PDG 2004 [18]. Equal 
production of neutral and charged B meson pairs in T(4S) decay is assumed. 
Reconstruction efficiencies are estimated by applying the same selection cri- 
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teria to 20,000 signal Monte Carlo events. The efficiencies for each final state 
are shown in Table 1. Here, the branching fraction of Kg tt+tt^ is 

not included in the efficiency calculation. In the fits to data, the width of the 
Xc2 is constrained to be 10% larger than the width of the Xci as expected from 
Monte Carlo simulation and the photon resolution. The effect of this assump- 
tion is included in the systematic errors. To calculate upper limits (U.L.), 
we add the statistical and systematic errors in quadrature, and assume that 
the combined error behaves as a Gaussian. We then follow the prescription of 
Feldman and Cousins [20]. The branching fractions are summarized in Table 
1. The first errors are statistical and the second errors are systematic. The 
statistical significance (E) of the signal in terms of the number of standard 

deviations is calculated as —2ln{CQ / Cmax) i where Cmax and Co denote the 
maximum likelihood with the nominal signal yield and with the signal yield 
fixed at zero, respectively. If a decay mode is reconstructed in two final states, 
e.g. B+ XciK*^{KIti+) and B+ XciK*+{K+7r^), the two branching frac- 
tion measurements are averaged, weighted by the statistical and systematic 
errors. Comparison between our measurement and those at BaBar [21,22] is 
also shown in Table 3. 



10 Systematic uncertainties 

Various sources of systematic uncertainty are summarized in Table 2. The un- 
certainty on the tracking efficiency is estimated to be 1.5% for e and n tracks, 
and 1.1% for /j, and K tracks. The pion tracks from Kg decay are from a 
displaced vertex and thus have a larger systematic error; we include 2% per 
track uncertainty for pion tracks from Kg decays. The errors on electron and 
muon identification are 1.6% and 2.2%, respectively. The kaon identification 
error is estimated to be 1.9%. 

Uncertainties on 7 and tt" reconstruction efficiency are 2% and 4%, respec- 
tively. The uncertainty in the Kg selection efficiency is checked by comparing 
yields for a sample of high momentum Kg — > 7r~^7r~ decays before and after 
applying the Kg selection criteria. The efficiency difference between data and 
Monte Carlo simulation is less than 1.0%. 

All modes have a systematic uncertainty due to Xci(2) lJ/'4^ and ^ 
t^l~ branching fractions; the uncertainties for these modes are taken from the 
PDG [18] and added in quadrature. 

The detection efficiency is found to have so small dependence as a function of 
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decay angles that the integrated interference term among amphtudes with dif- 
ferent quantum number of the i^'vr-system gets almost cancel. Therefore, this 
effect on the branching fraction is estimated to be less than 1.3%. However, 
the integrated amplitudes should be added coherently for the non-resonant 
Kt: background of the same quantum numbers. To estimate possible system- 
atic errors due to interference between the signal and background, we fit the 
i^Tyr-mass distribution to a function which assumes a "coherent" background 
contribution. We define this function as 



F{mK-.) = \^[N'sS{mK^) + ^J{N^)e'^Bg{mK-.)? + N^il - ae)S2(m^^)(l) 



is a relativistic P-wave Breit-Wigner function with tuk* — 892 MeV/c^, 
Bg(mKn) is the phase space background function, Ng and Nf, are the num- 
bers of events for signal and background, respectively, and Oc is the fraction 
of coherent background with a varied phase relative to the signal. The max- 
imum value of the coherent background fraction is taken from the helicity 
analysis to be Qc — 0.6. By varying the phase and relative fraction of coherent 
background within error in this fit, we find that the maximum difference in 
the signal yield due to the possible interference term is 15%. This difference 
is included as a systematic error and listed as "Inter." in Table 2. We assume 
that this difference is same for all K* channels; it is determined from the 
highest statistics B — > XciK* mode. 

In this study we have assumed that the background shape is described by 
the phase space function. We also fit the same distribution by using different 
background shapes. In one set of fits, a threshold function [17] is used to 
describe the background shape and in another set the function Bg{mK'n){ai -\- 
(I2PK + c^s-Pxci) describe the background shape; here BglmKn) is 

the phase space background function, ai, a2, 03 are the free parameters in 
the fit, Pk is the momentum of K in the K* rest frame and P^^^j is the 
Xci momentum in the B rest frame. The maximum difference in the signal 
yield obtained by using different background shapes is found to be 10.8% 
and is taken as an additional systematic error listed as background shape 
("Bkg. Shape") in Table 2. This difference is also determined from the highest 
statistics B — > Xci^* mode and assumed to be the same for all K* modes. 

The detection efficiency can depend on the various values of polarization am- 
plitudes for B Xci{2)K* decay modes. To check this effect, the exclusive 
Monte Carlo sample of B ^ XciK*{K~^'k~) decay with various values of the 
polarization amplitudes is generated. The maximum difference in detection 
efficiency is estimated to be 5.1% and treated as systematic error. This differ- 



where 



S{mKn) 



(2) 
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ence is assumed to be the same for all K* modes. 

The yield measurements are taken from fits where the mean and width of the 
signal shape are fixed. The corresponding uncertainty is estimated by varying 
the means and widths in turn by ±la, and repeating the fit. The percentage 
change in the signal yield is presented in Table 2. The difference between two 
background fitting models (i.e. using the J/ip inclusive background shapes 
from data or from the ip inclusive Monte Carlo) has also been included as a 
systematic error listed as background error ("Bkg. Error") in Table 2. 

For the Xci mode, the error on the efficiency correction due to the Xci width is 
based on the uncertainty of the fitted width in data. For the Xc2 mode, there are 
additional effects due to the assumptions in the Ht: a{xc2) '■ o-{xci) ~ 1.1 : 1 
(estimated in an alternative fit where the Xci and Xc2 widths float indepen- 
dently) and the uncertainty on the correction factor for the feed-across from 
Xci- The resulting error estimates are combined in quadrature, and listed as 
an efficiency correction ("Eff. Corr.") in Table 2. 

The uncertainty on the number of BB events, which is common to all modes, 
is 1%; this is negligible compared to the other systematic errors. 



11 Conclusion 

We have measured the branching fractions for the B Xci{2)K{K*) decay 
modes using 152 million BB events. Although inclusive Xc2 signals in B decay 
have been observed before [7], we do not find any significant signals for the 
exclusive two-body B — > Xc2K{K*) decay modes and set upper limits on 
the branching fractions at the 90% confidence level which are an order of 
magnitude smaller than those for B XciK{K*). The hclicity distribution 
for the B — > XciK* decay mode shows that the longitudinal polarization 
component is dominant. 
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Fig. 1. The (M^+^-^ — M^+£-) mass difference distribution with the selected Xci and 
Xc2 niass regions. 
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Fig. 2. The Mbc distribution for (a) B+ XciK+, (b) B+ Xc2K+, (c) 
XciKg and (d) B^ Xc2Kg decay modes, respectively, together with re- 
sults of the fit for the signal region and contributions due to J/ip inclusive and 
feed-across backgrounds. 
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Fig. 3. The i^7r-mass distribution for (a) XciK*^{K+tt-), (b) 

5° ^ Xcii^*°(i^g7rO), (c) B+ XciK*+{K+Tr°), and (d) B+ ^ Xcii^*+(i^g7r+) 
decay modes, respectively, together with results from the fit to K* signals and 
contributions of backgrounds due to J ftp inclusive, non-resonant and feed-across 
backgrounds. 
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Fig. 4. The i^vr-mass distribution for (a) 5° Xc2^*°(i^+7r-), (b) 

5° ^ Xc2if*°(ifg7r0), (c) 5+ ^ Xc2if*+(i^+7rO), and (d) S+ ^ Xc2-fC*+(if§7r+) 
decay modes, respectively, together with results from fit to K* signals and contri- 
butions of backgrounds due to J/ip inclusive, non-resonant and feed-across back- 
grounds. 
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Fig. 5. The efficiency-corrected cos Ok* distributions for Xci^*^i^~^'^~) decay 

mode fitted with F {cos Ok*) = A + 3Scos^ 9k* function for (a) background region 
(Mktt < 0.8 GeV/c^) and (b) signal region (0.821 GeV/c^ < Mr^t < 0.971 GeV/c^). 
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Table 1 

The summary of the estimated branching fractions for B Xci{2)^iK*) decay modes. The first ( 
quoted error is statistical (systematic). 



Mode 


Efficiencv 


Yield 


Br. Fraction (B) 


U.L.(90% C.L.) 


Sign. 




f(%) 




xlO-^ 


xlO-^ 


(S) 

v / 




25.4 


645.7 ± 26.9 


4.49 ± 0.19 ± 0.53 




37.4 




18.5 


127.0 ± 11.8 


3.51 ±0.33 ±0.45 




15.3 




17.4 


209.9 ± 22.7 


3.19 ± 0.35 ± 0.73 




12.1 




4.2 


5.0 ±4.8 


1.83 ± 1.76 ± 0.70 




1.1 








3.14 ±0.34 ±0.72 








11.4 


53.6 ± 11.1 


3.61 ± 0.75 ± 0.89 




6.9 


^ yciK*+(i^+7r°) 


7.0 


62.6 ± 13.2 


4.73 ± 0.99 ± 1.09 




6.0 








4.05 ± 0.59 ± 0.95 






-B+ Xc2K^ 


24.9 


14.8 ± 7.7 


0.16 ± 0.08 ± 0.02 


0.29 


2.1 




18.3 


2.6 ±2.1 


0.11 ±0.09 ±0.02 


0.26 


1.6 


5° ^ xc2i^*°(i^+7r-) 


16.5 


11.7 ±9.4 


0.29 ± 0.24 ± 0.08 




1.3 




4.6 


4.3 ±3.5 


2.24 ± 1.82 ±0.53 




1.6 








0.32 ±0.23 ±0.08 


0.71 




5+ ^ Xc2i^*+(i^°vr+) 


11.2 


5.2 ±5.4 


0.56 ±0.58 ±0.13 




1.1 


S+ ^ Xc2i^*+(i^+7r0) 


7.3 


2.5 ±7.2 


0.28 ± 0.82 ± 0.24 




0.3 








0.47 ±0.47 ±0.13 


1.27 





Tabic 2 

Summary of the multiplicative systematic errors in terms of percent (%) for all modes. Here the columns list the individual sys- 
-tematic error components: Trk., Lep. and Had. are those associated with tracking, lepton and hadron efficiencies; 7/7r° refers to 
7 and tt^ reconstruction; Kg refers to the Kg selection efficiency; B.F. refers to the errors on secondary branching fractions; Inter 
refers to the interference between K* signal and non-resonant Kn background; Bkg. Shape refers to different background shapes 
used in studying interference; Ang. refers to the error due to the uncertainty of the angular distributions; Fit. Para, lists sensiti- 
-vity to the fit parameters; Bkg. Errors refers to different fitting models for the J/^ inclusive background; and Eff. Corr. lists 



errors associated with the efficiency correction. 





Efficiency Uncertainty 


Jj.r . 


Inter . 


Bkg. 


Ang. 


r It 


Bkg. 




iotai 


Moaes 


Irk. 


Lep. 


Had. 


11 


Ks 






Shape 




TD „ f _\_-\ —.\ 
irQXS,.\±L(T ) 


Errors 


Corr. 






70 


07 
70 


07 
70 


07 
70 


07 
70 


07 
70 


07 
/o 


07 
70 


07 
/o 


Q7 
/o 


07 
70 


Q7 
70 


Q7 
70 


XciJ< 


3.7 


2.7 




2 




iU.O 










O.O 


O.O 


11. / 


^ XciK" 


D.D 






2 


1 


ID. 5 










0.9 


O.O 


12. y 


B — > Xcl-f^ (A^vr j 


5.2 


2.7 


1.9 


2 




iU.O 


15.0 


10. 


5.1 




Z.U 


O.O 


zz.y 


5° ^ Xci^*°(i^°7rO) 


6.6 


2.7 




6 


1 


10.5 


15.0 


10.8 


5.1 




30.0 


0.0 


38.3 


B+ ^ XciK*+{K^T,+) 


8.1 


2.7 


0.8 


2 


1 


10.5 


15.0 


10.8 


5.1 




7.0 


0.4 


24.6 


B+ ^ XciK*+{K+7r'') 


3.7 


2.7 


1.1 


6 




10.5 


15.0 


10.8 


5.1 




0.8 


0.3 


23.2 


B+ ^ Xc2K+ 


3.7 


2.7 




2 




8.5 










8.1 


2.0 


12.9 


^ Xc2K^ 


6.6 


2.7 




2 


1 


8.5 










7.7 


3.8 


14.2 


BO ^ ^^^K*^{K+^-) 


5.2 


2.7 


1.9 


2 




8.5 


15.0 


10.8 


5.1 




15.4 


1.7 


26.9 


5° ^ Xc2i^*°(i^°vr0) 


6.6 


2.7 




6 


1 


8.5 


15.0 


10.8 


5.1 




4.7 


2.3 


23.6 


B+ ^ Xc2i^*+(i^°7r+) 


8.1 


2.7 


0.8 


2 


1 


8.5 


15.0 


10.8 


5.1 


3.8 


1.9 


1.9 


23.2 


B+ ^ Xc2i^*+(i^+7r0) 


3.7 


2.7 


1.1 


6 




8.5 


15.0 


10.8 


5.1 


8.0 


84.0 


4.0 


87.4 



Table 3 

Comparison between this measurement and those at BaBar [21,22]. The first (second) quoted error is stat- 
-istical (systematic). 





Belle 


BaBar 




Br. Fraction xlO ^ 




Br. Fraction x 10 ^ 


TT T tono/ T ^ NX 1 n— 4 




4.4y ± u.iy ± u.oo 




K 70 -I- n 0(5 -I- n fid 
0. / y ± u.zd ± u.oo 






3.51 ± 0.33 ± 0.45 




4.53 ±0.41 ±0.51 






3.14 ±0.34 ±0.72 




3.27 ± 0.42 ± 0.64 




B+ ^ xcii^*+ 


4.05 ± 0.59 ± 0.95 




2.94 ± 0.95 ± 0.98 






0.16 ± 0.08 ± 0.02 


0.29 


0.09 ±0.10 ±0.11 


0.30 




0.11 ±0.09 ±0.02 


0.26 


0.21 ±0.11 ±0.13 


0.41 




0.32 ±0.23 ±0.08 


0.71 


0.14 ±0.11 ±0.14 


0.36 




0.47 ±0.47 ±0.13 


1.27 


-0.15 ±0.05 ±0.14 


0.12 



